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This paper develops an analytic theory for the study of some Poly a urns with random rules. The idea is to extend the 
isomorphism theorem in Flajolet et al. (2006), which connects deterministic balanced urns to a differential system for 
the generating function. The methodology is based upon adaptation of operators and use of a weighted probability 
generating function. Systems of differential equations are developed, and when they can be solved, they lead to 
characterization of the exact distributions underlying the urn evolution. We give a few illustrative examples. 
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1 Introduction 

A classical definition of the Polya urn scheme specifies it as an urn which contains balls of up to two 
colors (say black and white), and which is governed by a set of evolution rules. A step then consists in 
randomly picking a ball from the urn, placing it back, and depending on its color, adding or discarding a 
fixed number of black and white balls. 

Flajolet et al. (2006) introduces an analytic method to deal with a class of two-color Polya urns with 



deterministic addition rules. The method leads to a fundamental connection to differential equations via 
an isomorphism theorem. In the present manuscript, we extend the method to deal with Polya urns with 
random addition rules. One can think of schemes with random replacements as mixtures of deterministic 
schemes, with appropriate probabilities. 

We develop an analytic theory for the study of balanced tenable Polya urns with random rules (we 
shall make the terms "balanced" and "tenable" more precise in Section |2j. Balance and tenability are 
restrictions on the rules that admit an analytic treatment through generating functions. Our main result 
is the following: to each system of random rules satisfying the balance and tenability conditions, one 
can associate a differential system which describes the generating function counting the sequence of 
configurations of the urn scheme. A formal statement of this theorem appears in Section[3] 
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When the system can be solved, it leads to characterization of the exact distributions underlying the urn 
evolution. In Section [4j we give some examples where the generating function is explicit. Some of them 
include uniform and binomial discrete random variable in their rules. Other examples investigate some 
coupon collector variants. Finding the generating function and deducing the probability distributions is 
quite mechanical, owing to classical analytic combinatorics. 

2 Basic definitions 

In Polya's classical urn model, we have an urn containing balls of up to two different colors, say black 
and white. The system evolves with regards to particular evolution rules: at each step, add (possibly a 
negative number) black and white balls. These rules are specified by a 2 x 2 matrix: 

a,5eZ, /3, 7 eZ> . (1) 

The rows of the matrix are indexed with the colors black and white respectively. The columns of the 
matrix are indexed by the same colors: they are from left to right indexed by black and white. 

The asymptotics of this construct have been approached by traditional probabilistic methods in |Athreya| 
|and Karhn| ( |1968| ), [Smythe| ( |1996| l and |Janson| ( |2004| l. |Mahmoud| p008) presents some of these findings. 

In this manuscript we deal with a model where the constants a, (3, 7, and 5 are replaced by random 
variables. We attempt to provide an approach for the study of the urn's composition after a finite, and 
possibly small, number of ball draws. This is an equally important line of attack and can be viewed as 
more important in practice]®] when the number of draws is not sufficiently large to warrant approximation 
by asymptotics. We focus on exact distributions in this paper. 

2. 1 Balanced tenable urns with random entries 

We consider urns for which the dynamics of change are now represented by the replacement matrix 

{i *) • 

where A, B,C and T> are discrete random variables. If a random variable has a negative realization, it 
means we discard balls. 

Let B n and W n be respectively the number of black and white balls after n draws from the urn. We 
call the pair (B n , W n ) the configuration of the urn after n steps (draws). We start with a deterministic 
initial configuration (Bq, Wo) = (bo,Wo). At each step, a ball is uniformly drawn from the urn. We look 
at its color and we put it back in the urn: if the color is black, we add A black balls and B white balls; if 
the color is white, we add C black and T> white balls; the random variables are generated independently at 
each step. 

We now recall a few definitions, that serve to delineate the scope of our study. 

Definition 2.1 (Balance) An urn with the replacement matrix ([2]) is said to be balanced, if the sums across 
rows are equal, that is if A + B = C +T>, and these sums are always constant and equal to a fixed value 
9. The parameter 8 is called the balance of the urn. 

W For example, the Ehrenfest urn is a model for the exchange of gases. If the gas is being exchanged between an airconditioned 
room and the outside, the user is more interested in what happens within the next hour or so, not what ultimately happens at 
infinity. 
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All urns considered in this paper are balanced. As a consequence, the random variables are governed 
by the relations B = 9 — A, and C = 8 — T>. In addition, we will only consider urns with nonnegative 
balance, 8 > 0. The case 8 < refers to diminishing urn models. Questions and results on these models 
are completely different; see Hwang et al. ( 2007) for an analytic development. 



Since we allow the random variables A and V to take negative integer values, it is necessary to consider 
the concept of tenability. 

Definition 2.2 (Tenability) An urn scheme is said to be tenable, if it is always possible to apply a rule, 
i.e., if it never reaches a deadlocked configuration (that is, choosing with positive probability a rule that 
cannot be applied because there are not enough balls). 

To illustrate this notion, suppose we have an urn of black and white balls, with the replacement matrix 

-2 2 

2 -2 

This matrix, together with an even initial number of black balls and an even initial number of white balls 
form a tenable scheme. Any other initial condition makes this scheme untenable. For instance, if we start 
with two black balls and one white ball, the scheme gets stuck upon drawing a white ball. 

Remark 2.1 A tenable balanced urn scheme cannot have negative realizations for the nondiagonal en- 
tries B = 8 — A and C = 8 — T>, and they cannot have unbounded supports, either. Indeed, if B has a 
negative realization, say —k for positive k, we discard k white balls upon withdrawing a black ball. In 
view of the balance condition, we must add at the same step 8 + k black balls; we can apply this rule 
repeatedly until we take out all white balls, and the scheme comes to a halt (contradicting tenability). 
Therefore, B and C must have only nonnegative realizations. Furthermore, suppose the urn at some step 
has N black balls; if B has an unbounded support, it includes a value greater than 8 + N. Here A 
must realize a value less then —N. A draw of a black ball requires taking out more than N black balls 
(contradicting tenability). 

Example 2.1 This example illustrates a mixed model urn in the class we are considering. It has random 
entries in the replacement matrix and is balanced and tenable. Consider the Polya-Friedman urn scheme 
with the replacement matrix 



where B p is a Bev(p) random variable 



(ii) 



B p I -B p 
I -Bp B p 

Such a scheme has balance one, and alternates between 



(3) 



J Jj, and (J I 

that is between Polya-Eggenberger's urn ( with probability p) and Friedman's urn ( with probability 1 — p). 
The transition in the limit distribution of black balls when p goes from to 1 is quite unusual. For p = 
(Friedman's urn), a properly normalized number of black balls has the classical Gaussian distribution, 
andforp = 1 (Polya-Eggenberger's urn), a properly scaled number of black balls has a beta distribution. 
The scheme with replacement matrix (JiJ is an interesting example that warrants further investigation, and 
will be a special case of two illustrative examples (Subsections \4.2\ and \4.3\ . 

The notation Ber(p) stands for the Bernoulli random variable with parameter p, that assumes the value 1 with success probability 
p, and the value with failure probability 1 — p. 
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2.2 Definitions for the analytic methodology 

The first use of analytic combinatorics for the treatment of urn models is in |Flajolet et al.| ([2005 ), where 
some balanced urns are linked to a partial differential equation and elliptic functions. We now describe 



the analytic tools we use for our study. For a more complete (yet accessible) account see Morcrette ( 2012 1 



(§2). For a precise overview of analytic combinatorics, see the reference book Flajolet and Sedgewick 
( |2009] >. 

The total number of balls after n draws, denoted by s n := B n + W n , is deterministic because of the 
balance condition. Indeed, at each step, we add a constant number of balls (the balance 8), so 

s n = «o + On = 6 + w + On . (4) 

The main tool in this study is the weighted probability generating function 



Z 

Q(x, y, z):=J2 Yl s ° Sl ■ ■ ■ s "-! P( - Bn = b > W n= w ) xb V W — ■ (5) 



n=0 b.w>0 



The variable x counts black balls, y counts white balls, and z counts the number of draws. As b, w and n 
are governed by Q, one of the variables is redundant. So we can set y = 1 without loss of generality: it 
is sufficient to study a marginal distribution because it determines the joint probability distribution. 

Remark 2.2 The kernel sqSi . . . s n _i corresponds to the total number of paths of length n starting from 
the initial configuration. Indeed, at step we can pick a ball among the s balls in the urn; at step 1 we 
can pick a ball among the s\ balls in the urn, etc. 



Remark 2.3 The coefficient [x b y w z n /n\]Q(x,y, z) does not represent a combinatorial count; this 
quantity can be fractional. Nevertheless, it is closely connected to the notion of "history " described in 
\Flajolet et al.\ §2006 ). A history is nothing but a finite sequence of draws from the urn. One history of 
length n describes one possible evolution from step to step n. It can be viewed as a path in the evolution 
tree of the urn. In the context of deterministic rules, histories are combinatorial counts. However, in our 
case, we need more information to have a complete description of a sequence ofn draws. At each step, we 
can choose different rules depending on the realization of the random entries at that step. So, we are now 
dealing with weighted histories, that is tlie sequence ofn draws with a weight factor corresponding to the 
probability of applying the n different rules we use along the path from step to step n. If we denote by 
Qn,b,w the contribution of weighted histories beginning at step in the configuration (bo, Wq) and ending 
at step n in the configuration (b, w), we have 

oo n 
n=0 6,uj>0 

The kernel sqSi . . . s n _i is obtained by adding all weighted histories of length n, that is by setting x and 
y to 1: 

Q(l, 1, z) = (1 - Q Z y s °l e ■ nl [z n ]Q(l, 1, z) = SaSl . . . Sn _! = 9 n j' 9) . 

(in) \Ve use the operator [rr^ 1 . . . x J k k ] to extract the coefficient of x^ 1 . . . x J k k from a generating function of xi, . . . , x^. 
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In this way, we preserve the equivalence between the probabilistic model and the combinatorial aspects 
of weighted history. The balance condition is key to this equivalence: 

\x h y w z n ] Q(x, y, z) Q^hm 
F(B n = b,W n = w)= 1 ",J.*\ ' = ^ nAw • (7) 

[z n \ Q[l, 1, z) s ...s n -i 

3 An isomorphism theorem for urn schemes with random entries 

In this section, we state our main result through the following theorem which links the behavior of a 
tenable balanced urn to a differential system. The theorem always yields a differential system, and the 
difficulty then lies in extracting information on the generating function: indeed, the system obtained may 
not be always explicitly solvable by currently known techniques. 

Theorem 3.1 Given a balanced tenable urn with the replacement matrix 

A 6-A 
' — V V 

where A and T> are discrete random variables with distributions 

V(A = k)=n k , P(V = k)=r k , 
for —K < k < 9, (for some K > 0), the probability generating function of the urn is given by 

Q(x,y,z) = X b »(z)Y w °(z), 

where (bo, Wq) is the starting configuration, and the pair (X(t) , Y(t)) is the solution to the differential 
system 



*'(t) = E *kx k+1 {t)y e - k {t); 

k=~K 
9 

y'(t) = J2 r k x e - k (t)y k+1 (t), 



k=~K 

applied at t — 0, and any initial conditions x := a;(0) and y :— y(0), such as xy 0. 



Proof: The proof follows and generalizes that in Flajo let et al.| (|2006). The basic idea in this proof is that 



the monomial x l yi represents the configuration (_B n , W n ) — (z, j). The actions of the urn transform such 
a configuration into x l+k yi +e ~ k , if a black ball is drawn and if the nth realization of the replacement 
matrix is such that A = k, and 9 — A = — k (which occurs with probability Hk); there are i choices for 
such a black ball. Alternatively, the monomial is transformed into x l+e ~ k yi +k , if a white ball is drawn 
and if the nth realization of the replacement matrix is such that 6 — T> = 9 — k, and T> = k (which occurs 
with probability T k ); there are j choices for such a white ball. We represent the transformation by the 
operator 

e 

M:= J2 KkX k y e - k Q x +T k x - k y k Q v , 
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where 6 n = ud u is the pick and replace operator (lv) When this operator is applied to the configura- 
tion (i, j), it yields 



M(x i y j ) = J2 ilT k x i+k y i+e ~ k + j n x^"~ K y J 



r i+9-k J+k 



k=-K 



This operator is used to represent all possible configurations of the urn at step n + 1, given the configu- 
rations at step n. Indeed, if we write Q(x, y, z) := X^^Lo 1n(x, y)z n /n\ , as a function in one variable z, 
the operator A4 represents the transition 

q n+1 (x,y) = M (q n (x,y)) thus q n (x,y) = M n (q (x,y)) = M n (x b °y w °) . 

Let x(t) and y(t) be two functions that have Taylor series expansion near t = 0. Recall that the Taylor 
series expansion of the product x e (t + z) y m (t + z) is 



\t + z)y m (t + z) = J2 Q tn {At)v m {t)) 



z 
n ! 



We also have 



Q (x(t),y(t), z) = M " (^°(*) ^°(*)) 



n=0 



Z" 



Take £ = bo, and m = too. and these two expansions coincide, if the operators A4 and d/dt are the same. 
This is possible if, for all i, j > 0, 



Y in k x i+k (t)yi +9 - k (t)+jT k x i+9 - k (t)yi +k (t) = ix^y^t) x'{t) + j aty)^"" 1 y'(t) , 

k=-K 

and can happen by choosing 

( 

x'(t) = J2 KkX k+1 (t)y e - k (t); 

k=~K 
6 

y'(t) = Yl r k x e - k (t)y k+1 (t). 



k=-K 



Hence, if (X(t), Y(t)) is a solution of this system with initial conditions x := X(Q) and y := Y(0), we 
have Q (X(t),Y(t), z) = X b ° (t + z) Y w ° (t + z)0The statement follows by letting t = 0. □ 



(lv * See Subsection 1.6.2 p. 86 in Flajolet and Sedgewick 1 2009 1 for more details on the Symbolic method. 

(v ' A classic theorem of Cauchy and Kovalevskaya guarantees the existence of a solution, and the condition xy ^ guarantees an 
analytic expansion near t = 0; see, for example, jFolland ( 19951 (Chap.l Sect.D., p. 46-55). 



Analytic solvable balanced tenable urns with random entries 



7 



4 Examples of exactly solvable urns with random entries 

We now consider a variety of examples for which Theorem |3.1| admits an exact solution, and an exact 
probability distribution is obtained. In Subsection 4. 1 we treat a variation of the coupon-collector problem. 



In Subsections |4.2| and 4.3 we study urns with binomial random variables and uniform random variables. 
In Subsection |4.4| we study a variation of the coupon-collector problem that is modeled with three colors, 
demonstrating that the analytic methods can be extended to more colors. 



By probabilistic techniques, Janson (20041 and Smythe ( 1996) provide broad asymptotic urn theories 
covering some of the urns in the class considered here; our focus is on the exact distributions. 

4. 1 An urn for coupon collection with delay 

Collecting k coupons can be represented by a number of urn schemes. One standard urn model considers 
uncollected coupons to be balls of one color (say black) and collected coupons to be balls of another color 
(say white); initially all k coupons are not collected (all k balls are black). When a coupon is collected for 
the first time (a black ball is drawn from the urn), one now considers that coupon type as acquired, so the 
black ball is recolored white and deposited back in the urn. Collecting an already collected coupon type 
(drawing a white ball from the urn) results in no change (the white ball is returned to the urn). Thus, the 
standard coupon collection is represented by the replacement matrix 

-1 1 


Occasionally, the coupon collector may misplace or lose a collected coupon immediately after collecting 
it; the coupon still needs to be collected. That is, the drawing of a black ball may sometimes result in no 
change too, delaying the overall coupon collection. This can be represented by the replacement matrix 



13p I3p 











(8) 



where B p is a Ber(p) random variable. Under any nonempty starting conditions, this scheme is balanced 
and tenable. 

Proposition 4.1 For the urn the exact probability distribution of the number B n of uncollected 
coupons (black balls) after n draws is given by 



bo 



Sq-PJ 
so 



Proof: The urn scheme ([H) is amenable to the analytic method described, and leads to an exactly solvable 
system of differential equations. By Theorem |3.1| the associated differential system is 

f x'(t) = (l-p)x{t)+py{t); 

\ v'(t) = y(t), 

with solution 

X(t) = ye 1 + (x — y) e (1 " p)t , and Y(t) = ye 1 . 



8 



Basile Morcrette and Hosam M. Mahmoud 



Thus, 

Q(x, 1, z) = (e z + (x- l)e {1 - p)z ) bo e w ° z . 

In the course of coupon collection the total number of coupons (balls in the urn) does not change; that is, 
s„ = so, for all n > 0. Extracting coefficients, we get 

50 •••f"' 1 ¥(B n = b) = [x b z n ] e w ° z Y ( b °) (x - lV'e^-PW* x eS b °-^ z 
n - J 

/b °\ r_i v- b ( A e ((i-pW+6o-j)» 



0=b 

Since s = 6 + wo> we can write this as 



Si;)*- 1 * v„ 



Consequently, we find the probability expression. □ 

4.2 A binomial urn 

Consider an urn scheme with replacement matrix 



%6,p " — Xg p 
— ^0, P Xe, P 



(9) 



where Xg p is distributed like Bin(#,p) (vl) Under any nonempty starting conditions, this is a balanced 
tenable scheme of balance 6. For an expression of the probability distribution of black balls after n draws, 
it seems that a solution for general p is too difficult. We focus here on the tractable unbiased case p = J. 

Proposition 4.2 For the unbiased case p = ^ of the urn (|5J), the exact probability distribution of the 
number of black balls B n is 

1 / On \ 

¥(B n = &) = —(_ L for b G {b , b + l,...,b + On}. 

In other words, we have B n = bo + Bin (On, 1/2) . Asymptotics follow easily from the normal approxi- 
mation to the binomial distribution. Let Af(0, a 2 ) be a normal random variate with variance a 2 . We see 
that 



where the symbol — > stands for convergence in distribution. 



(vl ' The notation Bin(fc, p) stands for a binomially distributed random variable that counts the number of successes in k independent 
identically distributed trials with rate of success p per trial. 
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Remark 4.1 The replacement matrix Q is the special case where 9 = 1, and as mentioned in the in- 
troduction it alternates between Polya-Eggenberger' s urn (with probability p) and Friedman's urn (with 
probability 1 — p). So, one may think of the urn process as a mixture model that chooses between two 
given models at each step. It is well known that the number of white balls in a pure Polya-Eggenberger 
urn (scaled by n) converges to a beta distribution (see Pdlya ( 7937|J, wher eas the number of white balls 
in a pure Friedman's urn satisfies the central limit tendency (see Freedman {1965)): 



B,, 



In an unbiased mixture we get 7V(0, 1 /4) as limit, with centering and scaling similar to that in Friedman 's 
urn. The normal limit for an unbiased mixture has a larger variance than a pure Friedman 's urn, in view of 
occasional perturbation by the Polya-Eggenberger choice. Because Friedman's urn is "mean reverting," 
the entire mixed process is mean reverting, that is having tendency for average equilibrium around an 
even split. We thus see that tlie Friedman effect is stronger than the Polya-Eggenberger. 



Proof of Proposition |4.2| By Theorem 3.1 the associated differential system is 

Q 

X '(t) = £y(i- P )«-* (f)** + W-*(*); 

y'(t) = Y.P k ^-p) e - k ( 9 Xx e - k (t)y k+1 {t). 

k=0 ^ ' 

The unbiased case, p = 1/2, yields an explicit solution: 



X(t) = 



x + y\ 8 t ) 1/e 



and Y(t) = 



The probability generating function Q(x, y, z) = X(z) b ° Y(z) w " gives 



x + y \ 6 \ 1 / e 



(i-.(=±S)'«) 



Q(x,y,z) = ^2 ^2 s si...s n -iF(B n =b,W n =w)x b y v 



n=0 b,w>0 



x b °y w ° 



x + y \ 9 



6 0+ m 



Note that F(B n = b, W n = w) is for all values of b and w, except for w = On + so — b, when the 
probability may differ from zero. In other words, the joint probability distribution can be determined from 
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either marginal distribution. We compute 

Q(x, 1, z) 



s s 1 ...s n _ 1 ¥(B n = b,W n = w)x b 

n=0 b,w>0 
oo oo 

s si . . . s„_i 



n =b)x b 



n=0 6=0 





re! 



x 

'x + V 



(vii) 



Extraction of coefficients yield 

[a; b z"]Q(x,l,0) 



S ■ ■ • Sn-l 



S n = b) 

x+l 



[ x b z n ] x bo J2(-e( 



bo+w 

e 

n 



e„ S p ■ ■ ■ S ra -1 
n! 



0n 

— V 



0ra\ b s ■ • ■ s„-i 

X : . 



2 9n ^ V b 

6=0 v 

From this last expression we deduce the stated probability distribution for the number of black balls. □ 

4.3 A uniform urn 

Consider an urn scheme with replacement matrix 



u e e-u e 

9 — Mo Lia 



(10) 



where Ug is a uniformly distributed random variable on the set {0, 1, ... , 8}. Under any nonempty starting 
conditions, this is a balanced tenable scheme. 

Proposition 4.3 The exact probability distribution of black balls for the urn \10\ is given by 



F(B n = b) = T , e o n : h T^ , for b G {b , b + l,...,b + On), 



where T e ^ k := [x k ] (1 + x + . . . + x e ) n . 



The numbers Tg tn< f. are known in the classical literature (see Euler ( 1801 1) and have numerous com- 
binatorial interpretations. For instance, in Flajolet and Sedgewick (2009i (Subsection 1.15, page 45), 
T T —\ k n-k i s tne number of compositions of size n with k summands each at most r. In other words, 



< vil > Recall that for n, k > , 



-fc(-fc-l)...(-fc+ra-l) 
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T2, n ,k< called the trinomial coefficient, is the number of distinct ways in which k indistinguishable balls 
can be distributed over n distinguishable urns allowing at most two balls to fall in each urn. Several other 
interpretations related to counting strings, partitions and unlabeled trees of height 3 exist; see [Andrews] 
( fl99Ul l. 



Proof of Proposition |4.3| By Theorem 3.1 the associated differential system is 



x'(t) 

y'(t) 



k=0 



k=0 



This yields an explicit solution: 



X(t) = 



1 - 



e 



T G>v-'» 



1/0 



Y(t) 



1 - 



;=o i=o 
Thus, the probability generating function Q(x, 1, z) — X(z) b ° Y(z) w ° is 



1/0 



Q(x,l,z) 



oo oo 

sosi • ■ ■ s n -i 

n=0 6=0 



b) x b 



■III 



1 - 



tG>'W 



1=0 



Extraction of coefficients yields 



s ■ ■ ■ s n-l 

n! 



b) 



\x b z n ] x bo 



\x b ] x bo 



n=0 
1 



1=0 



(e + iy 



l + x- 



bp+Wg 



n Sq ■ ■ ■ S n -1 



□ 



We thus obtain the stated probability expression. 

4.4 An urn for a two-type coupon collection 

The method has obvious extensions for k > 2 colors. Consider, for instance, the balanced tenable three- 
color coupon collection urn scheme with the entries 



-1 B p I -B p 





(11) 



12 



Basile Morcrette and Hosam M. Mahmoud 



in which a collected coupon is randomly categorized to fall in one of two classes. Here again, B p is a 
Ber(p) random variable. The uncollected coupons correspond to black balls in the urn (replacement rules 
on the first row), and the white balls of the previous example are now ramified into red and green balls 
(second and third rows). The columns of the matrix are indexed by the same colors: they are from left to 
right indexed by black, red and green. Let B n , R n , G n be respectively the number of black, red and green 
balls in the urn after n draws. 

Proposition 4.4 For the urn the distribution of red balls R n is given by the following probabilities: 




Proof: This scheme has the underlying differential system 

x'{t) = py(t) + (l-p)h(t); 

y'(t) = v(t); 
h'(t) = h(t) . 

This system of differential equations has the solution 

X(t) = (py+(l-p)h)(e t -l) + x; 
Y(t) = ye'; 
H(t) = he 1 . 

Note that in this scheme the total number of balls after any number of draws remains the same at all times. 
The corresponding generating function is 

Q(x,y,h,z) = Y^ E 4nBn = b,R n = r,G n =g)x b y r W — . 

n=Q b,r,g>0 

The extended isomorphism theorem then gives 

Q(x,y,h,z)=X b °(t)Y r °(t)H^(t). 

One can extract univariate and joint distributions from this. For instance, 

s" 

-±V(R n = r) = [y r z n ]Q(l,y,l,z) 

= [V r z n ] (p (y - 1) (e z - 1) + e z ) bo y ro e roZ e 9oZ . 

Extracting coefficients (details omitted) we find the probability distribution for red balls. □ 
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5 Conclusion 



We presented a generalization of the isomorphism theorem for deterministic schemes. The generalization 
covers the class of balanced tenable two-color urn schemes with random entries. The result is particularly 
useful when the system of ordinary differential equations is amenable to a simple solution, from which 
exact probability distributions can be extracted in a fairly mechanical fashion. 



Our last example (see Proposition 4.4 1 shows it is straightforward to extend our theorem to additional 
colors: an urn with k colors will lead to a differential system of k equations. This theorem invites further 
work, using asymptotic tools from Flajolet and Sedgewick (2009), to obtain limit laws. For instance, the 



first example (see Proposition 4.2 a model mixing Polya and Friedman urn models, and Fig. [T]and|2]) 
raises various questions such as: Can we obtain an exact explicit expression for the generating function 
for general p e (0, 1)? What is the limit law, and what is the asymptotic behavior? How does the phase 
transition between a Gaussian regime and a beta regime take place? 






(a) p = 



(b) p = 0.4 



(c) p = 0.8 



(d) p = 1 



Fig. 1: Each graph represents a simulation of 100 histories of length 50 for the Polya-Friedman's urnJSl (see 



Example 2.1 1 with parameter p £ {0, 0.4, 0.8, 1}. There is a transition between a Gaussian behavior (Fig 1(a) I and 



i beta behavior (Fig |l(d)[ . 




Fig. 2: Normalized distribution of the number of black balls after 200 draws (B200) for the Polya-Friedman's urn ([3} 
with different values of p G {0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1}. Here, the starting configuration is one 
black and one white ball. We observe a transition between a Gaussian distribution and a uniform distribution (which 
is a special case of beta distribution, when bo = wq = 1). 
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